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Abstract. In this paper, we focus on the thin film equation with lower order 
"backwards" diffusion which can describe, for example, the evolution of thin 
viscous films in the presence of gravity and thermo-capillary effects, or the thin 
film equation with a "porous media cutoff' of van der Waals forces. We treat 
in detail the equation 

u t + {u n {u xxx + vu m ~ n u x - Au M ~ n u x )} x = 0, 

where v = ±1, n > 0, M > m, and A > 0. Global existence of weak non- 
negative solutions is proven when m — n > —2 and A > or v = — 1, and 
when —2 < m — n<2, A = 0, v = 1. From the weak solutions, we get strong 
entropy solutions under the additional constraint that m — n > —3/2 \{ v = 1. 
A local energy estimate is obtained when 2 < n < 3 under some additional 
restrictions. Finite speed of propagation is proven when m > n/2, for the case 
of "strong slippage," < n < 2, when v = 1 based on local entropy estimates, 
and for the case of "weak slippage," 2 < n < 3, when v = ±1 based on local 
entropy and energy estimates. 

Keywords: thin film equation, backwards diffusion, higher order parabolic equations, 
degenerate parabolic equations, finite speed of propagation 
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1. Introduction 



The thin film equation [44] 

u t + {u n (u xxx )} x = 0, n > 0, (1.1) 

with n = 3 models the dynamics of thin viscous films with no slip boundary condi- 
tions, and with n = 1 it models Hclc-Shaw flow. Often (1.1) needs to be augmented 
with various lower order terms in order to take into account the presence of addi- 
tional physical effects, and certain such equations shall be considered here. See [42] 
for a survey and review. Some of the physical systems which are accommodated by 
the analysis in this paper include: 

(i) the evolution of thin viscous films in the presence of gravity and thermo-capillary 
effects 

u t + {u n (u xxx + u m - n u x - Au M - n u x )} x = 0, (1.2) 
where m, n, M, A are constants such that A > 0, n > 0, m < M, as well as the 
more accurate variant of (1.2) given by 

u t + {u n (u xxx + h'(u)u x )} x = 0, (1.3) 

where h'(u) = —vG + M ( 14 ^ 2 „)2 , G,Bi,B% are positive constants, V = ±1, and 
where V = +1(— 1) represents stabilizing (destabilizing) gravitational forces [42, 52, 
43]. Equation (1.3) with A = 0, 

u t + {u n (u xxx + u m - n u x )} x = 0, 

models thin films with thermo-capillary effects but without gravitational effects 
when m — n — — 1 [21], or with destabilizing gravitational effects but without capil- 
lary effects when m — n = 0. The value of n in (1.2), (1.3) reflects the assumptions 
on the slip conditions at the interface of the thin film with the underlying substrate, 
with n = 3 modeling no slip and < n < 3 modeling various types of slip. 

(m) the equation 

u t + {u n (u xxx + h'{u)u x )} x = 0, (1.4) 
describes (i) the evolution of a thin viscous film in the presence of attractive polar 
forces if h(u) = —b\e~ u l h ^ where bi, hi are positive constants, or (ii) the evolution of 
a thin viscous film in the presence of attractive van der Waals forces if h(u) = Bu~ b , 
where B < is a (negative) Hamaker constant and b is a positive constant. More 
generally, if lim„^o^'(w) > (< 0), h{u) is said to represent limiting attractive 
(repulsive) forces. Equations of the form (1.2), (1.4) can represent a combination 
of attractive and repulsive forces, and the limiting power m = lfrn^o uh"(u)/h'(u) 
can assume a range of both positive and negative values in modeling various forces 
such as polar forces, van der Waals forces, as well as "porous media cutoff" of van 
der Waals forces. See [38, 42, 11, 31, 16]. 
(Hi) the Hocherman-Rosenau equation [29] 

u t + {f{u)u xxx + g{u)u x )} x = 0, (1.5) 

was proposed as a generalization of the cylindrical Kuramoto-Sivashinsky equation, 
which models low Reynolds number two phase cylindrical flows. Equation (1.5) 
with f(u), g(u) > has been used as a prototype equation for studying the relative 
strength of the second and fourth order terms in determining criteria for blow up 
[29, 14]. Setting f(u) = u n and g(u) = u m or g = h'(u) with h'(u) > in (1.5) 
yields (1.4) with "limiting attractive forces." 
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To understand these various models, we shall focus on the equation 

u t + {u n (u xxx + vu m - n u x - Au M - n u x )} x = 0, (1.6) 

where v = ±1, and where n, m, M, A are constants satisfying A > 0, n > 0, 
M > to, with some further restrictions to be imposed in the sequel. All of the 
examples outlined above may be written in the form (1.6), except for (1.3), (1.4 i) 
in which 

h'(u) = u(l + Bu)- 2 or h'{u) = ^e- u/b \ (1.7) 
respectively. Noting that 

u 1 

0< (TTb^ < W (L8) 
<|l e -«/fe < ^ ; (L9) 

02 2 

for u > 0, we see that the examples in (1.7) have upper bounds of the form 
h'(u) = vu m ~ n - Au M ~ n with v = 1, A = 0, and m - n = -1 and to - n = 
respectively, which facilitate the analysis of (1.3), (1.4 i). Our treatment of (1.6) 
can be generalized to encompass (1.7) as well; some remarks in this direction are 
included within the text. 

In terms of the physical systems being modeled, it is reasonable to assume zero 
contact angle, u x — 0, and zero flux, u n {u xxx + vu m ~ n u x — Au M ~ n u x ) = 0, along 
the external boundaries of the domain, and thus that 

u x = u xxx = if u > and x = ±o. (1-10) 

While these boundary conditions are adopted here, other boundary conditions are 
possible to consider. The zero contact angle condition reflects the physical assump- 
tion that the viscous liquid film completely wets the underlying substrate; non-zero 
contact angle conditions, while physical, have proven up to now to be difficult to 
implement in dynamical problem formulations [41]. A simple alternative to (1.10) 
is to impose periodic boundary conditions, as we shall amplify. 

The term " {u n {vu m ~~ n u x )} x " with v = +1 in (1.6), is often referred to as a 
" backwards diffusion" term, since if one considers dynamics dominated by this term 
alone 

u t + {u n (u m - n u x )} x = 0, 

and one linearizes about a uniform positive state, then the resultant dynamics is 
given by the backwards (ill-posed) diffusion equation. Similarly, if v = —1, the 
term '{it" '{vu m ~ n l u x )} x " in (1.6) is often referred to as "forward diffusion," for 
obvious reasons [40]. We shall often refer to equation (1.6) with v = +1 as the 
"unstable case" and to equation (1.6) with v = — 1 as the "stable case," since in 
the context of thin films, (1.6) with v = +1 models limiting attractive forces which 
are destabilizing ("long wavelength unstable" in the terminology of [13]) and (1.6) 
with v = — 1 models limiting repulsive forces which are stabilizing [42] . 
The degenerate Cahn-Hilliard equation 

u t + {u(l - u)[- Inu + ln(l - u) + au + e 2 u xx ] x } x = 0, (1-H) 

whose solutions satisfy < u < 1 [23, 40], can also be said to be of the form (1.6) 
with v = — l,n = l,m = for u near or 1. Equations similar to (1.4), (1.6) also 
arise in modeling structure formation and the dynamics of biofilms [33] , as well as 
in modeling the dislocation density in plasticity theory [26]. 



Equation (1.6) in the presence of forward or stabilizing diffusion 

u t + {u n (u xxx - u m - n u x )} x = 0, (1.12) 

has been studied somewhat more thoroughly than the backwards or unstable variant 

u t + {u n (u xxx + u m - n u x )} x = 0. (1.13) 

This is perhaps not surprising, as the behavior of the thin film equation changes 
qualitatively 1( ?ss upon c idding a stabilizing term as opposed to adding a destabilizing 
term, and hence the analytical tools used in studying (1.1) may be more readily 
adapted to its analysis. 

Since (1.1) is designed to model the evolution of this viscous films, nonnegative 
initial date should yield nonnegative solutions. Compactly supported initial data 
whose support which may spread or shrink are also physically relevant to consider. 
A pioneering step in the analysis of (1.1) was the proof of the existence of weak 
nonnegative solutions for 1 < n < 3 by Bernis & Friedman [3] for nonnegative H 1 
initial data, which relied on energy and entropy estimates. The analysis in [3] was 
extended by Bernis [4] to the interval < n < 3. In [8], existence (nonexistence) 
was demonstrated of compactly supported spreading source type solutions to (1.1) 
for < n < 3 (n > 3), and at the contact line, the solutions were seen to behave 
like (xo — x) 2 for < n < 3/2 and like (xo — x) 3 /™ for 3/2 < n < 3. By developing 
refined entropy estimates, the existence of strong (C 1 for a.e. t > 0) nonnegative 
solutions was demonstrated by Beretta, Bertsch & Dal Passo [1] and by Bcrtozzi & 
Pugh [12], and the solutions were seen to possess the regularity of the source type 
solutions at the contact line. Positivity properties were seen to depend strongly on 
the value of n [3, 1], with "touchdown" being possible for small values of n [10, 1]. 
Solutions were proven to become positive in finite time, and to converge to their 
mean as t — > oo [1, 12]. The only nonconstant steady states for (1.1) are of the 
form (x — a) + (b — x)+ for a < b, which do not possess the regularity of the strong 
solutions [1, 12]. Finite speed of propagation of the support of strong solutions for 
< n < 3 was proven in [5, 6]. Some of the details of the various methodologies 
will be amplified further as we present our results. Since these basic studies, the 
analysis has been developed considerably to encompass, for example, the Cauchy 
problem for measured valued initial data [18], waiting time phenomena [17], higher 
dimensional studies of (1.1) [19, 9], as well as various numerical schemes [28, 27]. 

With regard to (1.12), existence of nonnegative distributional solutions was 
proven in [11] for < n and < m < 1. These solutions were shown to approach 
their mean as t — > oo, becoming positive in finite time. Formal asymptotics and 
numerics were used to suggest the existence of advancing fronts for n > 3, as well as 
for < n < 3. Existence of nonnegative solutions for f2 C R N , N = 1, 2, 3 was then 
demonstrated for two special cases of (1.12); namely, for the model for defects in 
plasticity theory mentioned earlier [26], and for (1.11), the Cahn-Hilliard equation 
[23]. In [20], an existence theory for (1.12) was presented which encompassed the 
cases n > 1/8, m > — 1, N < 3 (with n < 4 if N = 3), with the case < n < 1/8 
being treatable by adjusting the definition of the solution. The existence results in 
[11, 20] build on the results in [3, 1, 12], and in [20] it builds on [26] as well. In 
particular, use is made of augmented entropy estimates [1, 12], which indicate that 
the contact angle is zero for the solutions obtained if < n < 3 and n—m > — 2. The 
constraint n — m > — 2 allows the energy to be bounded from below for arbitrary 
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nonnegative initial data Uo 6 H 1 . Regularity at the contact line was confirmed 
via formal asymptotic arguments in [11]. In [20], finite speed of propagation is 
demonstrated for m > 0, 1/8 < n < 2, and estimates are obtained for the speed of 
propagation; infinite speed of propagation is proven form < 0, —2 < m — n < —3/2. 

In [2], self-similar spreading source type solutions were shown to exist for (1.12) 
when m — n = 2, if < n < 3, and not to exist if n > 3. The asymptotics 
of these solutions at the contact line was seen to match that of the source type 
solutions found for (1.1) in [8]. In terms of steady states, in [34] nonconstant 
positive periodic steady states and zero contact angle steady states are shown not 
to exist, though touchdown steady states with non-zero contact angle are shown to 
exist if m — n > —2. 

In [28, 27], existence of positive solutions from nonnegative initial data was 
demonstrated via construction of a numerical scheme, for (1.6) with v = — 1 and 
m — n < —2. The inclusion of the term —Au M in (1.6) can be expected to enhance 
the rate of decay of the solution to its mean, and should not have any major effect 
on contact angle and propagation properties. In our analysis in the present paper, 
we closed gaps in the parameter range for (1.12) and (1.6) with v = — 1, both in 
terms of existence and the finite speed of propagation property. Moreover, local 
energy and entropy estimates are obtained. 

With regard to (1.13), less attention has been paid to the various existence, 
regularity, and finite speed of propagation properties for (1.13) than for (1.12), 
but there has been longstanding interest in the qualitative predictions of (1.13), in 
particular in regard to the possibilities of rupture [44, 49] and blow up [29]. There 
has been considerable interest in various self-similar and steady state solutions for 
(1.13), and their stability [47]. 

In terms of existence, regularity, and finite speed of propagation properties for 
(1.13), the results up to now may be summarized as follows. In [13] it is demon- 
strated that if < n < 3 and n < m < n + 2, then there exist globally bounded 
nonnegative weak solutions to (1.13) which possess the finite speed of propagation 
(FSP) property, and have the same regularity at the contact line as was found in 
[8] for (1.1). We remark that formal asymptotics developed in [11] for (1.12) are 
equally valid for (1.13). The methodology in [13] relies on [3, 1, 12, 11], incorporat- 
ing a "disjoining pressure" potential arising from the lower order terms [31] into the 
energy used in the basic energy estimate. Numerical evidence is given that blow 
up may occur if m > n + 2. The case n = 1, m>3 = n + 2is considered in [14], 
and it is demonstrated that for compactly supported nonnegative initial data whose 
initial energy is negative, there exists a solution which blows up in finite time. For 
the critical case, n = 1, m = 3, this implies blow up if the initial mass is sufficiently 
large [51]. It was conjectured in [13, 14] that m — n > constitutes a necessary 
condition for well posed dynamics. The existence proof in [27] mentioned earlier in 
the context of (1.12), includes the case v = 1, i.e. (1.13) with n > and m < n — 1, 
for positive initial data. The inclusion of a term of the form —Au M in (1.6) when 
v = 1 should eliminate the possibility of blow up, and should not have any major 
effect on contact angle and propagation properties. 

The set of steady states and self-similar solutions is much richer for (1.13) than 
for (1.12). For (1.13), positive periodic steady states exist, as do compactly sup- 
ported "touchdown" steady state solutions with zero as well as with nonzero contact 
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angles [34] . A study of steady states and stability based energetic criteria was un- 
dertaken in [36] for m — n e [1, 2). It has been demonstrated [34, 35, 37, 47] that 
steady states with zero contact angle exist for all < n < 3, n < m, which are 
stable if m < n + 2, < n < 2, marginally stable if m < n + 2, 2 < n < 3, and 
unstable otherwise. In [15] for v = 1, n = 3, m = —1, steady states are seen to con- 
verge to a (5-distribution in the limit in which repulsive forces are neglected; these 
results should perhaps be compared with the non-single- valued profiles seen in [38] 
which result when the term u xx is replaced by the mean curvature, a correction 
which becomes important in the singular limit. In terms of self-similar solutions, 
both spreading and blow up self-similar solutions are possible. Rupture self-similar 
solutions have also been observed and studied [53, 50]. If m = n + 2, there ex- 
ist spreading self-similar solutions if < n < 3, and none if n > 3 [2]. Blow up 
self-similar solutions were considered in [46] , with existence being demonstrated for 
< n < 3/2 and non-existence for n > 3/2. For a discussion of the stability of 
self-similar solutions, see [51, 47]. 

The focus of the present paper, however, is not on blow up, steady state solu- 
tions, and self-similar behavior, but rather on conditions that guarantee existence, 
regularity, and finite speed of propagation, completing and enhancing what was 
previously know, in order to construct a framework in which understand, for ex- 
ample, the transition between rupture, positivity and touchdown properties, global 
existence, and blow up. Roughly speaking, with regard to existence for v = — 1, 
< n < 3, there had been a gap for n — 2 < m < — 1; and now we have existence 
of both weak and strong solutions for n — 2 < m < —1, as well as for n — 2 = m 
for a constrained set of initial conditions. For v = 1, < n < 3, previous existence 
results had required m > n, and now we have existence of weak solutions in the 
interval n — 2 < m < n, as well as strong energy /entropy solutions in the subin- 
terval n — 3/2 < m < n. Note also that for v = ±1, < n < 3, in the interval 
n — 2 < m < n — 1, we do not require that the initial data be strictly positive [15]. 
In terms of FSP, for v = —1, < n < 3, there had been a gap at < n < 1/8 
which has now been filled for < m < n + 2, as well as a gap at 2 < n < 3, which 
has now been filled for n/2 < m < n. For v = 1, < n < 3, there had been a gap 
at < m < n, which has now been filled for n/2 < m < n. 

What physics pertains to the interval v = — 1, < n < 3, n — 2 < m < — lor 
v = 1, < n < 3, n — 2<m<n? The thin film equation with v = 1, < n < 3, 
— 2<m — n<0 can reflect a "porous media" cutoff of attractive (or repulsive if 
v = — 1) van der Waals forces [42], thin films under the influence of thermocapillary 
effects with n = 2or3, m — n = —1, v = 1, [21] (with other values of n, < n < 3, 
also being possible to consider), or a restricted Hocherman-Rosenau equation [29]. 
If one takes v = l,n = 2or3,m< —1, which can model, for example, the effects of 
a Lennard- Jones potential or a thin film on a layered solid substrate in the limit in 
which the limiting repulsive forces are neglected, the dynamics are known to lead 
to rupture [49, 53, 50]. 

As a first step in this direction (enhancing existence), the existence of weak non- 
negative solutions (see Definition 1) is demonstrated in §2. This is accomplished 
by means of the basic energy estimate [3, 1, 12] 




T 



(1.14) 
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where il = (—a, a), a G (0, oo) is arbitrary, C depends only on the problem pa- 
rameters and the initial conditions, and < T < oo. The estimate (1.14) holds for 
< n, — 2 < to — n, to < M, with the additional constraint that to — n<2ifz/=l 
and A = 0. Additionally, an entropy estimate [3, 1, 12] is obtained, using Gron- 
wall's inequality for regularized solutions. The use here of Gronwall's inequality, 
which explicitly depends on the regularization parameter, differs somewhat from 
elsewhere [28], allowing us to control the lower order, possibly singular, forcing 
terms. These estimates, together with mass conservation, 

/ u(x,t)dx= / u(x, 0)dx, (1-15) 
Jn Jn 

imply global bounds from which existence of weak nonnegative solutions based on 
uniform Holder continuity and continuation arguments. See Bernis & Friedman [3], 
Giacomelli [24]. 

To obtain the existence of a strong (C : (fi) for a.e. t > 0) solution, a local 
entropy estimate is derived in §3, following [1]. For v = 1, the additional constraint 
to — n > — | is imposed and the local entropy estimate may be written as 

1 



f CV+ Q (x, T)dx + A [ C 4 u a+M - 1 u 2 x dxdt+ 
a{a + 1) J n J Qt 

[ ( 4 u a+n - 2 ~< +1 (u r )l x dxdt+ [ ( 4 u a+n - 3 u x dxdt < 
J p Jq t 

c 2 f (\C x \ 4 + \CG x \ 2 )u n+a+1 dxdt + c 3 [ \(( 3 G) x \u a+m+1 dxdt+ 

JQt JQt 

c 4 I ( 4 u a+2m - n+1 dxdt+—^ f ( 4 u a+1 dx, (1.16) 

Jq t a(a + 1) J n 

where Q T = fl x (0, T), P = ft x (0, oo) \ {u = or t = 0}, and which holds for 
certain a G (max{— 2to + n — 1, —to — 1}, 2 — n) \ {0, —1} and for 7 satisfying 
(3.8). For v = — 1, a similar estimate is obtained without additional restrictions. 
For both v = ±1, the local entropy estimate implies the global entropy estimate 
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f u a+n - 3 u x dxdt + c e [ u a+n - x u xx dxdt<ciT, < T, (1.17) 

JQt JQt 

which holds also for (1.1), and which implies strong solutions and certain positiv- 
ity properties [1, 12]. We also present a refinement of Theorem 3.1 from [1] (see 
Corollary 3.4 in §3), which clarifies the set of (5 for which C 1 ([— a, a]) regularity 
for almost every t > is implied for u 1 /^(-, t) based on the local entropy estimates 
and the properties of the initial data. As in [25], a local energy estimate is derived 
for 2 < n < 3, under the additional restrictions that to > (2n — 2)/3if2<n<5/2 
and to > n - 3/2 if 5/2 < n < 3. 

In §4 and §5, we investigate conditions for the finite speed of propagation prop- 
erty (FSP) for (1.1). Two different techniques have been developed for studying the 
FSP property which both rely on energy and entropy estimates. The first method 
was developed by F. Bernis [5, 6] for the standard thin-film equation, (1.1). For 
< n < 2, Bernis in [5] introduced the weighted integral entropy function 

c T / 2 

E T {r):= / (r- |a;|) 4 (u ii± f ±i ) dxdt, < r < r , < T, (1.18) 

J0 J-r ^ ' xx 
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for < r < a, which corresponds to the third term on the left hand side of (1.16) 
with ((x) = (r— |x|)+. By taking suppuo C fi\{|a;| < r } and deriving a nonlinear 
ordinary fourth order differential inequality for Et{t), it was possible to demon- 
strate the FSP property. For 2 < n < 3, Bernis demonstrated the FSP property in 
[6] based on a similar ordinary differential inequality for an appropriately defined 
weighted integral energy function. 

An alternative method for studying propagation properties for various thin- 
film like equations was proposed for the case < n < 2 in [32] (sec also the 
references and comments in [30]). This method relies on functional rather than on 
differential inequalities for various entropy and energy functions. The emphasis in 
this technique is on cut off functions rather than on weight functions, which lead to 
functional dependence on parameterized subdomains of Cl. The main feature of the 
inequalities which are derived is that the minimal power of u(x,t) in the terms on 
the left hand side of the inequality is strictly less than the minimal power of u(x, t) 
which appears in the terms on the right hand side. Such inequalities in conjunction 
with the functional Stampacchia Lemma (see Lemma 6.1) and its generalization for 
the system of functional inequalities (see Lemma 6.2) can be used to demonstrate 
the FSP property in various contexts. In the present paper, we shall rely primarily 
on this second method. 

This latter method (the Fl-method) has been used [30] to obtain sharp estimates 
on the speed of propagation of the support of solutions to equation (1.1) with 
arbitrary u e Li(Q), ft \ suppu ^ 0. In [20], using this method for the thin film 
equation with forward diffusion (1.12) with source type initial data, m crit = n + 2 
was shown to be a transitional value between the long and short time asymptotics 
of the Barenblatt porous media behavior and the predicted behavior of source 
type solutions for the standard thin film equation [8]. Using the FI- method, sharp 
sufficient conditions were obtained in [17] relating the flatness at the edge of the 
support of m to the phenomenon of a waiting time for propagation of the support 
of solutions. Moreover for equation (1.1) with an additional nonlinear absorption 
term, in [45] a sharp sufficient condition on the flatness of uo guaranteeing the onset 
of shrinkage of the support was obtained, as well as estimates from below on the 
shrinkage rate. In future work, we hope to undertake similar analyzes for (1.6) with 
v=\. 

For (1.6) with "forward diffusion" (with A = 0), conditions for FSP as well as 
sharp estimates for the speed of propagation were obtained in [20] for the " strong 
slippage" case (0 < n < 2) only. The results which we obtain here are for (1.6) with 
v = — 1, with 2 < n < 3, as well as for the much more delicate case of backward 
diffusion, v = 1, where the lower order diffusion term "encourages" the destruction 
of the FSP property at all values of to. Our analysis makes use of some ideas from 
[25]. The proof given in §4 is for the "strong slippage case" in which < n < 2, 
and requires that m> njl'xiv = \. It relies on the local entropy estimate from §3 
for a positive and the Stampacchia Lemma for systems. We demonstrate that if 
suppu C {x < 0}, then there exists a continuous function, s(t) satisfying s(0) = 0, 
and a positive time To such that suppu(-, t) C [—a, s(t)], s(t) < a V£ < To, and 
s(T ) = a. The proof in §5 is for the "weak slippage case" in which 2 < n < 3, 
and requires that to > n/2. It is based on combining local entropy estimates for 
— 1 < a < with the local energy estimates from §3, and again makes use of the 
Stampacchia Lemma for systems. We conjecture that the restriction to > n/2 for 
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v = 1 is sharp or close to sharp. Note that in the stable case v = — 1, the value 
m = ^ is critical in the context of asymptotics near the edge of traveling wave 
solutions with constant speed of propagation (see [11]). Since the asymptotics in 
[11] are valid also for the unstable case v = 1, we expect that FSP-property fails 
in the unstable case for some m near to ^. Thus we suspect expect that for v = 1, 
the restriction m > n/2 is sharp or close to sharp. 

As to further questions and future directions, a question of interest in the present 
context is to understand exactly how bad the backwards diffusion can become while 
maintaining something of the smoothing properties of the thin film equation with 
n > 0. Thus one should like to identify transitional values in terms of rupture, 
positivity, touchdown, and infinite as opposed to finite spreading rates. Finally, 
one should like to extend all aspects of the analysis to higher dimensions. 

The outline of the paper is as follows. The existence of weak non-negative solu- 
tions is proven in §2. Existence of strong energy-entropy solutions is demonstrated 
in §3. Finite speed of propagation is proven in §4 for the case of weak slippage and 
in §5 for the case of strong slippage. 



2. Weak solutions 

In this section, we follow Bernis & Friedman [3], relying on local parabolic 
regularity theory [22, 24] to obtain global existence of weak solutions, defined below. 
We regularize the initial data and use Gronwall's inequality in the context of the 
entropy estimate. This allows us to avoid regularizing the lower order terms [13, 14], 
and allows us to widen the range of validity of the results. 

Notation. Let = (—a, a) where a e (0, oo) is arbitrary, Q t — £1 x (0, t), 
< t < oo, and set P t = Q t \ {u = or t = 0}, Q = Qoo, P = -Poo- 
Let us consider the problem 

' u t + {u n {u xxx + vu m - n u x - Au M - n u x )) x = 0, {x, t) e Q T , 
(P) I u x (±a, t) = u xxx (±a, t) = when u(±a, t)^0, t G (0, T), 
, u(x, 0) = uo(x), x G fi, 
where v = ±1 and < T < oo. We shall assume the initial conditions to satisfy 

u e H^ri), u > o, u ^ o. (2.1) 

While we shall look for solutions on a finite interval, we can always consider the 
parallel Cauchy and periodic problems obtained by extending the initial conditions 
via periodicity and reflection. This will allow us, for example, to directly implement 
the generalized Bernis inequalities for nonnegative periodic functions, [25, Lemma 
B.l] in obtaining local energy estimates in §3. See also the remarks in [1, 14]. 

Definition 1. A function u g C ' 1 / 2 ' 1 /^ x [0, oo)) ni°°([0, oo); H 1 ^)) is said 
to be a weak solution of (P) if: 

(a) ueC 4 -\P), u>0, 

(b) u x (x, t) = u xxx (x, t) = when u(x, t) ^ 0, for (x, t) e <9S1 x (0, oo), 

(c) J = u n {u xxx + vu m - n u x - Au M - n u x ) e L 2 (P), 



(d) for all <p G Lip(f2 x (0, oo)) with compact support, u satisfies: 

I u(f> t dxdt+ / u n (u xxx + vu m ~ n u x - Au M ~ n u x )(t> x dxdt = Q, (2.2) 
Jq Jp 

(e) u(x, 0) = u (x) for x G ti. 

Given this definition, we formulate 
Theorem 1. For the following range of parameter values: 

(i) v = -1, < n, < A, n - 2 < m < M, 

(ii) u=l,0<n, 0<A, n-2<m<M, 
(Hi) u = 1, 0< n, = A, n-2<m<n + 2, 

there exists a solution to (P) in the sense of Definition 1 for w satisfying (2.1). 

We shall find a weak solution to (P) as the limit of a subsequence of smooth 
positive solutions to a regularized problem, with regularized initial conditions, uo e . 
We shall require that for e > and for some A G (0, 1), 9 G (0, 2/5], uo e satisfies 

u 0e eC^ x (Cl), < (±a) = <'(±a) = 0, u + e < u 0e < u a + 1, 

u eQ — > mo in H ((—a, a)) as e — > 0. 

Following [3], we set 

f*( s ) = ill, 4 ; (2-4) 

4 

when < s « 1, / e (s) «^if0<n<4 and f e (s) w |s|" if n > 4, which shall allow 
us to guarantee positivity of the approximants u e for t > 0. Here and in the section 
which follows, c i} di denote positive constants that are independent of e, and Ci(t) 
denotes a positive increasing function defined on (0, oo) that is independent of e; 
q, di, Ci(t) may depend on 17, uq, and the problem parameters, and their value 
may change from line to line. 

Proof. Let us define the approximating problem (P e ) 

' u t + {Mu)(u xxx + uu m - n u x - Au M - n u x )} x = 0, (x, t) G Q T , 

(P e W u x (±a,t)=u xxx (±a,t)=0, t e (0, T), 

u(x, 0) = u 0e (x), x G £1. 

Problem (P £ ) possess a unique maximal positive solution, u £ , such that u e G 
C 4,a, a/4 x [ Q) T£ > 0; gee p2, Theorem 6.3, p 302] as well as the remark 
following the proof given there. That the associated Cauchy problem with period- 
ically reflected initial data maintains the periodicity and reflection properties, can 
be seen by translating and reflecting the solution, then invoking uniqueness of the 
solutions to the Cauchy problem. 

Testing (P e ) with <f> = 1 and recalling (2.1), (2.3), it follows that 

< W e {t) = ubT < u + 1, t G (0, T), (2.5) 

where v := J^ v - The equality u^(t) — uq^ in (2.5) expresses mass conserva- 

tion. Note that (2.5) also holds for solutions to similarly defined approximating 
problems for (1.3), (1.4 i). 
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Setting 



h(s) = < 



m — n+ 



r-^. m,M^{n-l}, 



v In s — A 



M-n+l ■ 



m = n — 1, 



. ^ — -^--Alns, M = n-1, 
and testing (P e ) with — u exx — h(u e ), we obtain the energy estimate [3, 23, 31] 



/ 

Jn 



u e l - H(u e ) 



dx 



fe{u e )(u exxx + h'(u e )u ex ) 2 dx dt = 



rl 



dx, 



(2.6) 



where 



H(s) 



As" 



(m— n+2)(m-n+l) {M -n+2)(M -n+l) ' 



i/(s In S — s) 



As" 



(m — n+2) (m—n+1) 



(M-n+2)(M-n+l) ' 

A(s In s — s), 



m,M ^ {n - 1}, 
m = n — 1, 
M = n — 1 . 



Let f = 1 and —1 < m — n < 2. Then noting that 

- 1 < sins - s < s 2 , 0< s, (2.7) 
recalling the Gagliardo-Nirenberg inequality 

IMIlp(q) < cilKII^^IMI^) +c 2 ||u|| L i(n), 1 < P < 4, 
and that the approximating solutions are positive and satisfy (2.5), it follows that 

[ H(u e )dx <- f u t ldx + c 3 (f u e dx) = 7 / u e 2 dx + c 5 . (2.8) 
Jo 4 V n \Jq / 4 Jq 

Suppose the restrictions on the parameters stated in Theorem 1 hold and that 
— 2 < m — n < — 1 if v = 1. Then recalling (2.7), and noting that for Ai, A 2 > 0, 

< a < 8, 



and that by Holder's inequality, if < a < 1, then 



we obtain that 



/* \u\ a dx<a [ \u\dx+{l-a)\Q\, 
Jn Jn 

/ i?(u e ) da; < c 6 (u + 1) = c 7 . 



(2.9) 



From (2.6), (2.8), (2.9), we may now conclude that 

- ( u 2 x {x,t)dx+( f e (u t )(u txxx +vu t m - n u ex -Au e M - n u ex ) 2 dxdt < c 8 . (2.10) 
4 Jo JQ* 



Remark 2.1. TTie constraint that m — n < 2 w/ien 1/ = 1, A = has been imposed 
to guarantee global existence. If is = 1, A = 0, andm — n = 2, (2.10) remains valid 
ifuo is sufficiently small, and global existence is again implied. See [13]. 
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From (2.5), (2.10), we obtain that 

IKIU°°(o,t ; jri(n)) - C8 ' \\f} /2 ( u txx + h(u e )) x \\mQ t ) <c 9 . (2.11) 

The following argument from [3] is by now standard. Noting that u t = —J x , (2.11) 
implies that 

IMIl^O,*;^- 1 ^))' \\ J \\L 2 (Q t ) < C W , (2.12) 

and the estimates (2.11), (2.12) can be seen to imply the uniform Holder estimate 

||Ue|| C o, 1/2, i/8(q;) - Cl1 ' (2.13) 

See [3] for details. 

Remark 2.2. It is also possible to implement the above discussion when u 11 is 
replaced by f(u) in (F), for f(u) 6 C(0, oo) — > M + , and to work with [3] 

f i n f(u)u 4 
ef(u) +u 4 

under suitable assumptions on f(u). 

We now demonstrate roughly as in [3, 24] that u e > 4er > in ft x [0, r £ ], where 
cr = er(e). On [0, r e ), recalling (2.13) we know that < u e (x, t) < A, for some A 
which is independent of e. Multiplying the equation in (P e ) by G' e (u e ) [39, 3], where 



Ge(s) = - I g e (r)dr, g e (s) = - f 

J S J S 



dr 



Mr)' 
and integrating, 

f G e (u e (x,t))dx + [ [u e 2 xx - {vu e m - n - Au e M - n )u e 2 x ]dxdt= [ G t (u 0e )dx. 
Jn JQ t Jo. 

(2.14) 

If m — n ^ — 1, then integrating the term — vu™ n u 2 x by parts, 

G e (u e (x, t))dx + u 2 xx dxdt + A u t M ~ n u 2 x dxdt 
JQt JQt 

<\ I u 2 x dxdt+- 1 ——r^f u e 2m - 2n+2 dxdt+ [ G e (u 0e )dx, (2.15) 

z JQt Z{m-n+L) j Qt J n 

recalling (2.11) 

G e (u e (x, t)) dx + f —u 2 xx dx dt < c\2 I u e ~ 2 dxdt + f G e (uo e )dx. (2.16) 
JQt 2 JQt Jn 

2m — 2n + 2 

If m — n= —1, the term (lnu £ ) 2 replaces ^-n+iy 1 m (2-15). Then noting that 
ln 2 (s) < ci 3 s -2 , < s < A < oo, 

the estimate (2.16) again follows. 
Since by (2.4), 

it now follows easily that 

/ G e (u e (x, t)) dx + I \-u 2 xx dxdt <— / G e (u e ) dx dt + c 15 t + j G e (u 0t )dx. 
Jn jQ t 2 e Jq t J n 
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es- 4 <-^— 0<s, (2.17) 



Using (2.3), (2.17), we find that f n G e (uo e ) dx < c\q. Now we may use Gronwall's 
inequality to conclude that 

/ G e (u £ (x, t))dx < D £ (t) < oo, t G [0, t £ ). (2.18) 

where for all < e <C 1, D e (t) is a positive increasing function of t defined on 
[0, oo). As in [3], (2.18) can be seen to imply positivity. 

The solution, u £ (x, t), may now be extended to exist globally, as in [3, 24]. Select 
f £ {s) G C 2 (K) such that f £ (s) = f £ (s) for s > 2a, and f £ (s) > f £ (a) for all set. 
Thus u e (x, t) also constitutes a weak solution of 

Uet + {fe(u £ )(u £xxx + h'(u £ )u £x )} x = 0, 
satisfying the same initial and boundary conditions as before. For x £ € [—a, a), set 

v £ (x, t)= u e (£, t)d4- / f £ (u £ (x £ , 9)){u exxx + h'(u £ )u £x }(x £ , 9) dO. 

Jx e JO 

The regularity and positivity of u e (x, t) imply that v e (x, t) is well defined in D = 
Cl x (0, r e ) and satisfies 

| V £t + MMX'tyiVexxxx + VU^^Vc^- Au £ M - n V <ixx } =0, 

I w e (± a, t) = v exx (± a, t) = 0. 

Using parabolic regularity results for v e , enhanced regularity may be obtained for 
u e and hence for the w e -dependent coefficients in (2.19). Returning again to (2.19), 
additional regularity is obtained for v e , which allows us to conclude that u € (x, r e ) G 
C^ x (0,). Therefore the solution may be continued, in contradiction to the assumed 
maximality of the solution. 

Having demonstrated the global existence of positive approximants, u e , existence 
of a sequence {u £k } converging uniformly to a solution of (P) on Q T , for all < T < 
oo, as £fe — > 0, is now implied by Arzela-Ascoli and the uniform Holder estimates, 
(2.13), as in [3]. □ 

Remark 2.3. It can be readily verified that the results of Theorem 1 remain valid 
when m — n = —2 and v = — 1, if the initial conditions satisfy the additional 
constraint that J n H(uo) dx < oo. 

Remark 2.4. The existence of weak solutions for (1.3), (Hi), withh'(u) replacing 
vu m-n _ A u M-n j n j} e jl n in on ^ can fr e concluded for initial data satisfying (2.1), 
(2.3), by verifying that (2.8) holds for (1.3), that (2.9) holds for (Hi), and that 
(2.16) holds for both (1.3) and (H%), then arguing as above. 



3. Strong entropy-energy solutions 

To get strong entropy-energy solutions, we derive local entropy estimates [1], 
which give us strong solutions [1, 12], then derive a local energy estimate [25]. In 
obtaining the local entropy estimates and strong solutions, we follow [1] closely. 
Throughout this section the parameters will be assumed to satisfy the conditions 
in Theorem 1. Moreover, in referring to solutions of (P) and (P e ), we shall assume 
that u satisfies (2.1) and u 0e satisfies (2.3). Some further restrictions shall be 
introduced in the sequel. 
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Before deriving the entropy estimates, we present a lemma, which constitutes 
a refinement of Theorem 3.1 in [1], and which is useful for concluding regularity 
results from entropy estimates. 

Lemma 3.1. Let u{x 1 1) be a weak solution of (P) obtained as the limit of a subse- 
quence of solutions u e (x, t) of (P £ ). Suppose that for some a e — n, 2 — n), there 
exist constants a, c 2 , and S > 0, which do not depend on e, such that 

f <+"-^ +1 (^)L^dt< Cl , (3.1) 



and 



for all 7 satisfying 



JQt 



u« +n - 3 u e 4 x dxdt<c 2 , (3.2) 



1 + n + a 1 + n + a 
3 < 7 < 3 1- 8, (3.3) 

then w 1//3 (-, t) e ^([—a, a]) for all (3 E (0, „ + ^ +1 ) /or afonosf ever?/ t > 0. 

Proof. For any < /3 < K+ ^ +1 , we may choose 7 satisfying (3.3) such that < 
/?7 < 1. Setting q = 4 — ( 1 +^+ a ) anc i arguing as in the proof of [1, Lemma 3.1], it 
follows from (3.1), (3. 2) that for almost every t > there exists a C\{t) < 00 such 
that 

if u(y, t) = for some y E [—a, a], then 

Ku 7 )^ 4 ^ 79 ^, < Ci(£)|a; - y\ {q - 1)/q for ir e [-a, a}. (3.4) 

From (3.4), we find by integrating that for almost every t > 0, there exists a 
C2(i) < 00 such that 

if u(y, t) = for some y E [—a, a], then 

u(x,t) < C 2 {t)\x - y|^+i+T for a; G [-a, a]. (3.5) 

Since < fS'f < 1, we may combine (3.5) and (3.4) to obtain that for almost 
every t > 0, there exists a C^{t) < 00 such that 

if u(y, t) = for some y G [—a, a], then for x £ [—a, a], 

Ku^Or, i)| < C 3 (t)|x-»|^|a:-y|^ < <7 3 (i)|a: - j/^, (3.6) 

where n = i - 1 > and C 3 (t) < 00. □ 

From Lemma 3.1, two simple but useful corollaries follow. 

Corollary 3.2. [1] Let u(x,t) be a weak solution of (P) obtained as the limit of a 
subsequence of solutions u e (x, t) of (P £ ). Suppose that for some a 6 (i — n, 2 — n), 
i/iere exisi constants ci, C2, and <5 > 0, which do not depend on e, such that for all 
7 satisfying (3.3), the estimates (3.1) and (3.2) hold. Then u(-, t) 6 C x ([— a, a]) /or 
almost every t > 0. 

Proof. If a e (§ - n, 2 - n), then e (1, 2). Hence 1 e (0, ^_). □ 

Remark 3.3. Lfu(, t) E C 1 ([— a, a]) /or almost every t > 0, iftera u(x, £) is said to 
&e a strong solution in the sense of Bernis [5] . 
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Corollary 3.4. Let u(x,t) be a weak solution of (P) obtained as the limit of a 
subsequence of solutions u e (x, t) of (P £ ). Let "J denote a nonempty subset of (\ — 
n, 2 — n). If for all a € ^, £/iere exist constants c\, C2, and <5 > 0, w/wc/i do not 
depend on e, such that for all 7 satisfying (3.3), the estimates (3.1) and (3.2) hold, 
then w 1 / /3 (-, t) £ C 1 ([— a, a}) for all j3 E (0, „ +in f ^ +1 ) f or almost every t > 0. 

Proof. The result is an immediate consequence of Lemma 3.1. □ 

We now derive our primary entropy estimates. 
Let 

C € C 4 ([— a, a]) with support in (—a, a) and £ > 0, (3.7) 
or C = 1, and let [1, 12] 



es »+n-3 s a+l 



G e (3)= , ; " , ; - + 



(a + n - 4) (a + n - 3) a(a + l)' 

where a G (1/2-n, 2-n)\{0, -1}. Using C 4 <^(u e ) to test (P £ ) on Q T =flx(0, T), 
< T < 00, and treating the terms which also appear in the classical thin film 
equation as they were treated in [1], we obtain that for any 7 satisfying 

f + V(t-2)(l-2t) _ c t+l + y/(t-2)(l-2t) ^ 
3 3 
where t = a + n, there exist positive constants, C3, C4, which do not depend on e, 
such that 

/ C 4 G e K(x, T))d.x+ 

c 3 / C 4 u e Q+ " _27+1 «)L^*+ / ( 4 u £ a+n - 3 U<ix dxdt < 

JQt JQt 

[ C 4 G e (u 0e )dx + c 4 f (\( x \ 4 + \(( xx \ 2 )u? +a+1 dxdt + I, (3.9) 
Jn Jq t 



where 



/:= - / ( 4 ge(u e ){f e (u e )(vu e m - n u ex -Au e M - n u ex )} x dxdt, g e (u e ) := G' e (u e ). 
Integrating / by parts, 

1= f ( 4 9UUe)fe(Ue)(vur- n U e l-Au e M - n U e 2 x )dxdt+ 
JQt 

[ 4( 3 C x ge(u e )Mu e )(vu e m - n u ex - Au e M - n u ex )dxdt := I a +I b . 

JQt 

The term I a may be written as 

I a = [ C 4 (^ e " +m_1 - Au e a+M ~ x )u 2 x dxdt. (3.10) 

JQt 

For v = —1, both terms in (3.10) are non-positive. For v = +1, we estimate 
I a <6 [ ( 4 u f a+n - 3 u f idxdt + c 5 (S) [ C A u e a+2m - n+1 dxdt- 

JQ T JQt 

A f C i u t a+M - 1 u e 2 x dxdt, (3.11) 

JQt 



JQt 
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where S > is arbitrary. 

With regard to h, integration by parts gives that 



4(C 3 C), 



dx dt. 



' g e ( S )f e ( S )[vs m - n -As M - n ]ds 
As noted in [1], \g t (u £ )f £ (u £ )\ < c 6 u™ +a . Thus, recalling (2.11) and that M > m, 

h<c 7 [ \{t\ x ) x \u a + m+1 dxdt. (3.12) 
Jq t 

If v = — 1, we may combine the estimates on I a and lb to obtain 

/ ( 4 G t (u £ (x,T))dx + A [ ( 4 u* +M ~ 1 u e 2 x dxdt+ 
Jn Jq t 

c 3 f C 4 Ue a+n ~ 2l+1 {WJ)l x dxdt+ [ ( 4 u t a+n - 3 u £x dxdt + 

JQt JQt 



^ u *+™-i Ue 2 dxdt < f ( 4 G t (u 0e )dx + II, (3.13) 
Jn 



where 



n = c J (\c x \ 4 + \cc xx \ 2 K +a+1 dxdt + c 7 [ \(( 3 c x ) x K +m+1 dxdt. 

JQt JQt 

Similarly, if v = +1, the estimates yield 

f ( 4 G t (u t (x,T))dx + A[ ( 4 u e a+M - 1 u e 2 x dxdt+ 
Jn Jq t 

c 8 f { 4 u e a+n - 2 -< +1 (uZ) 2 xx dxdt+ [ C 4 u t a+n - z u 4 x dxdt < 

Jq t Jqt 

+ c 5 f (* Ut a+2m - n+1 dx dt + [ ( 4 G t (u 0e ) dx + II. (3.14) 
Jq t Jn 

To obtain bounds from (3.13), (3.14), we impose certain conditions on a and on 
the initial data. 

Remark 3.5. Suppose that uq satisfies (2.1) and < n < 3. Defining 



a = 



-1, 



< n < §, 
| < n < 3, 



(3.15) 



we see that a* + 1 > 0, a* £ [| — n, 2 — n), and J n ( 4 Uq +1 dx < +oo /or all 
a e (a*, 2 - n) and ( 6 C 4 ([-a, a]). 

In consideration of the above remark, we define 

Definition 2. Suppose that uo satisfies (2.1) and ( G C 4 ([— a, a]). T/ien we define 
ao(C) = inf a swc/i i/iat a > ^ — n and 

j n C 4 «o +1 < +°° -1, 
J a C 4 | lnuo| dx < +oo i/ a = — 1. 
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Remark 3.5 and the definition of Qo(C) imply that if < n < 3, then 

i — n < a (() < a* < 2 — n. 



(3.16) 



With regard to the stable case, the theorem below follows essentially as in [1, 20]. 

Theorem 2. (The stable case.) Suppose that v = —1, < A, < n, — 2 < m—n, 
with m < M if A > 0, and with J n uo dx < oo if —2 = in — n. 

i) Let 13 e (0, /3 ) where j3 = max j n+ 3 Q+1 , m+< j, 0+1 j and a = a (C = !)• T^e™ 
m 1 / ,3 (-, t) G C 1 ([— a, a]) /or almost every t > 0. 

m) Lei £ satisfy (3.7). Then, for any a <G (max{«o(C)j — m — 1}; 2 — n)/{0, — 1} 
and for any 7 satisfying (3.8), 

( 4 u 1+a (x, T)dx+ 
f C, i u a+n - 2 ~< +1 {u^)l x dxdt+ f ( 4 u a+n - 3 u 4 x dxdt 

JP JQt 

Q 4 u a+m ~ 1 u 2 x dxdt + A f Cu^^uldxdt < 

JQt 

+ \(( xx \ 2 )u n + a+1 dxdt + c 3 f \(( 3 ( x ) x \u a + m+1 dxdt+ 

Jqt 



a(a + 1) J n 



ci 



+ 



C2 



/ (IC 



^— / C A u a+1 dx. (3.17) 



Proof. Part i) follows by setting £ = 1 in (3.13), implementing Lemma 3.1, and 
noting the Holder regularity implied by the boundedness of Jq t (u 5 ) 2 dxdt. 
Part ii) follows easily from (3.13) by letting e — > and noting that a + m + 1 and 
a + n + 1 are positive in the indicated parameter range. □ 

Remark 3.6. If m < and m — n < — |, then the Infinite Speed of Propagation 
Property holds for the solutions discussed in Theorem 2 due to the boundedness of 



Jo {u £L± ^ ±1 ) 2 x dxdt. See [20, Corollary 2.1]. 



Theorem 3. (The unstable case.) Let v = 1, < A, < n, — § < m — n, 

m - n < 2 if A = 0, and m < M if < A. 

i) Let oto = ao(C = 1)> a i = max{ao, —2m + n — 1}, and fto = n+ 3 l+1 ■ Then 
u 1/,3 (-, t) G C 4 ([-a, a]), for all (3 6 (0, /3 ) /or aZmosi eijery i > 0. 

ii) For any £ satisfying (3.7), let a 2 = max{ao(£), —2m + n — 1, — m — 1}. Then, 
for any a E (a 2 , 2 — n)/{0, —1} and for any 7 satisfying (3.8), 

, 1 - / ( 4 u 1+a (x, T)dx + A { ( 4 u a+M - 1 u 2 x dxdt+ 
a(a + l)J n J Qt 

f ( 4 u a+n - 2 ~* +1 (ut) 2 xx dxdt+ f ( 4 u a+n - 3 u 4 x dxdt < 

\ 4 + \(( xx \ 2 )u n+a+1 dxdt + c 3 f \(( 3 ( x ) x \u a + m+1 dxdt+ 

JQt 



Cl 



C2 



/ (ICx 
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/?o> 3 3 "1 (3.19) 
1 d < n < 3. 



c 4 / C 4 " Q+2m_n+1 dxdt + ^i [ ( 4 u Q a+1 dx. (3.18) 

7q t a(« + 1) J a 

Proof. To prove part i), note that the condition — | < m — n implies that ol\ G 
[| — n, 2 — n), and that ol\ > ao, ai + 2m — n + 1 > 0, then implement Lemma 3.1 
with C = 1- Part m) follows by noting that «2 G [\^ n i 2 — "), a > ao, a+m+l > 0, 
and a + 2m - n + 1 > 0, then letting e ->• in (3.14). □ 

Remark 3.7. TTie results given in Theorem 3 also hold for the exceptional cases 
(1-3), (1.4 i), with A — and with m — n assuming the values m — n = —1 and 
to — n = 0, respectively. This can be easily demonstrated by following the arguments 
above, once one notices that estimates (3.11), (3.12) also hold for (1.3), (Hi) 
when the value of to — n is taken as —1 or 0, respectively, by utilizing the bounds 
(1.8), (1.9). 

Remark 3.8. When < n < 3, Corollary 3.2 and (3.16) imply that under the 
assumptions in Theorems 2 and 3, the weak solutions obtained as the limit of a 
subsequence of {u e (x,t)}, the solutions to (P e ), are in fact strong solutions in the 
sense of Bernis [5] . In particular this implies that solutions with compact support 
have zero contact angle at the edge of the support. 

Using (3.15)-(3.16), one can check that in the context of Theorems 2 and 3, 

2, 0<n<§, 
2 < 71 < 

These bounds correspond to the bounds which were obtained in [1] for the thin film 
equation (1.1). The regularity indicated in (3.19) also corresponds to the contact 
line regularity of the self-similar source type solutions which were demonstrated in 
[2] to exist for (1.6) when < n < 3, v = ±1, A = 0, and to = n + 2. In this 
sense the regularity obtained is sharp, though additional regularity is provided by 
the terms Jq t (u^^ )1 dx, Jq t A(u^t )^dx, for some parameter values. 

Remark 3.9. For < n < 3, with regard to the terms in (2.2) which do not 
depend on to or M , since the regularity obtained is the same as regularity which 
was obtained for (1.1) in [1], these terms can be interpreted in the same sense as 
in [12], which is distributional for | < n < 3 and is in the weaker sense of (2.2) 
for < n < |. With regard to the terms in (2.2) which depend on to and M, these 
terms may be readily seen to be interpretable in a suitable distributional sense if 
to > — 1 and in the weaker sense of (2.2) for m < — 1. Details of the case v = —1, 
< to < 1 are discussed in [11]. 

In the case of "strong slippage," in which < n < 2, the local entropy estimates 
provided by Theorems 2, 3 can be used to prove the finite speed propagation prop- 
erty for the strong solutions obtained there when to > ^, see §4. In the case of 
"weak slippage" in which 

2 < n < 3, (3.20) 
the finite speed propagation property for standard thin film equation (1.1) and Uq G 
was proved by F. Bernis [6] using local energy estimates, which were derived 
there. The proof of the energy estimates is based on certain integral inequalities for 
smooth functions which are positive on some interval, see [7] . Our proof of the finite 
speed of propagation property for (1.6) for to > ^ in the case of "weak slippage" 

18 



(3.20), is based on combined usage of local entropy estimates, which were derived in 
Theorems 2 and 3, and local energy estimates which are derived below in Theorem 
4. The derivation of the local energy estimates makes use of a generalization of 
the Bernis inequalities mentioned earlier, for nonnegative periodic functions, see 
Lemma 6.3. Somewhat similar methodologies were employed in [30, 25] where local 
energy estimates were also derived, and which also relied on combined usage of 
local entropy and energy estimates. Notably, the combination of local energy and 
energy estimates used here and in [25] relies on using local entropy estimates with 
— 1 < a < 0; previous approaches relied on local entropy estimates with a > in 
proving the finite speed of propagation property [1, 12, 20]. 



Theorem 4. Suppose that v = ±1, 2 < n < 3, < A, and that m < M if A > 0, 
and m < n + 2 if A — and v = 1. Suppose moreover if v = —1, that m— |ra > —1, 
and if v = 1, that m — |n >— |«/2<n<| and m — n > — | i/| < n < 3. 
Under these assumptions, the strong solutions obtained in Theorems 2, 3 satisfy the 
following local energy estimate 

f ( 6 \u x (x,T)\ 2 dx + d 1 f ( 6 ((u n i 1 f x + (u^ 1 )l x + (u^) 2 xxx )dxdt+ 
Jn Jq t 

di I ( 6 u n u 2 xxx dxdt < / C 6 \u 0x (x)\ 2 dx+ 

jQ T n{u>o} Jn 

d 2 I u^i^dxdt+^t^dxdt- 

Jqt 

(yu m u x - Au M u x )(u x ( 6 ) xx dx dt, (3.21) 



/ 



/Q T n{M>o} 

where ((x) is an arbitrary nonnegative function from C 4 ([— a, a]). 

Proof. In order to derive the local energy estimate, we test the equation in the 
approximating problem, (P e ), with — (£ 6 u ex ) x , and easily deduce that 



u ex (x,T)\ 2 dx + / ( 6 f e (u e )\u exxx \ 2 dx dt 

JQt 

Ic 



[ £ 

in 2 

fe(u e )u exxx [2u exx (( 6 ) x + ^(C 6 )^] dxdt- 



Qt 

[ f e (u e )(vuT~ n u ex - Au™~ n u ex )(u ex ( 6 ) xx dxdt+ 

JQt 



2 



/ ^\u 0ex \ 2 dx. (3.22) 
Jn 



Assuming (3.20) and noting (3.15), (3.16), it is easy to check that max{a*, — m— 
1} < in Theorem 2 and that max{a*, —2m + n — 1, — m — 1} < in Theorem 3. 
Hence if n satisfies (3.20), then for arbitrary ( satisfying (3.7) or £ = 1, a may be 
chosen to be fixed and to satisfy 

- 1 < a < 0, (3.23) 

in addition to satisfying the constraints indicated in Part ii) of either theorem. 
For a satisfying (3.23), we have due to (1.15) 
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J u e {x,T) 1+a dx <\Q\- a (^J u t (x,T)dx^ 

\n\~ a (^J u 0e dx^ <d 3 . (3.24) 

Setting £ = 1 in (3.13) and employing (2. 3), (3. 24), we obtain the following inequal- 
ity when v = — 1, 



C3 



f u^ +n - 2 -' +1 (uJ)l x dxdt+ [ u? +n - 3 u e 4 x dx 

JQt JQt 



dt 



+ 



f ( u a+m-i +Au? +M - 1 )u e ldxdt < [ G t (u 0t ) dx + d 4 < d 5 . (3.25) 

JQt JO 

Similarly, setting £ = 1 in (3.14), employing (2. 3), (3. 24), and noting that the as- 
sumptions on a imply that a > a 2 > — 2m+n— 1, we obtain the following inequality 
when v = 1, 



C'8 



:dt 



+ 



f u* +n - 2 ^ +1 {ulf xx dxdt+ f u* +n - 3 u ex dx t 
JQt JQt 

[ Au a + M ~ x u? x dxdt< [ G e (u 0e )dx + c 5 [ u« +2m - n+1 dxdt + d e <d 7 . 

JQt JO JQt 

(3.26) 

Next we pass to the limit e — > in (3.22). First, setting £ = 1 in (3.22) yields 
that for any 6 > 0, 

2- 1 / |u e;c (x,T)| 2 dx+ / f e {u e )\u exxx \ 2 dxdt < 2- 1 / |u 0e;c | 2 ^+ 
Jo jq t Jo 

(l+^|sup Ue | M - m ) / / e K)<-"|w ex |K xxx |rfa ; ^<2- 1 / \u 0ex \ 2 dx+ 



sf f e (u e )\u exxx \ 2 dxdt+ [ u« +n - 3 u ex dxdt+c(8) [ u 4 e m - 2n - (a+n)+3 dxdt. 

JQt JQt JQt 

(3.27) 

Suppose that 

Am - 2n - (a + n) + 3 > 0, (3.28) 

then setting S — 1/2 in (3.27) and using the estimates (3.25), (3.26), we deduce 
that 

/ \u<: x (x,T)\ 2 dx + / f e (u € )\u exxx \ 2 dxdt < d 8 . (3.29) 

JO JQt 

Remark 3.10. In the context of Theorem 2, when n satisfies (3.20), m > n—2 > 0, 
hence max{a*, — m — 1} = a* = — 1. Thus for arbitrary initial data satisfying 
(2.1), (2. 3), (3.28) is satisfied for some admissible a, if 

3 

m — —n > —1, 

which is stronger than the previous constraint, m — n > — 2. In the context of 
Theorem 3, when n satisfies (3.20), m > n — 2 > 0; hence max{a*, — m — 1} = 
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a* — — 1, but a there must also satisfy a > —2m + n — 1. it is easy to c/iecA; fftai 
(3.28) ZioZds /or some admissible a, i/ and cm/y if m, n satisfy the condition 

3m-2n> -2. (3.30) 

Recalling the constraint m — n > — | m Theorem 3, it is easy to check that (3.30) 
constitutes an additional constraint if n < |. 

Using the estimates (3.29), (3.25), (3.26), it is easy to check that the integrals 
on the right-hand side of (3.22) are uniformly bounded with respect to e if (3.28) 
is satisfied. For arbitrary n > 0, u e — > u strongly in the space C 4:1 ({u > n}). 
Therefore passage to the limit e — > in all of the integrals in (3.22) over the 
domain {u > n} is straightforward. As to integrals over the domain {u < i]}, we 
have, for example, by virtue of (3.28), 



< 



/ fe{Ue)u exxx U ex U e m n ( 6 dx dt 

JQ T n{u<ri} 

1/2 , . v 1/4 



/ f e (u e )\u exxx \ 2 dxdt) ([ u« +n - 3 u tx dxdt) 

JQ T n{u<r)} / \JQ T n{M<Jj} / 

„ \ 1/4 _ 2 

/ u 4m-2n-( a +n)+3 ^ j < ^ 4,n-2n-(a+n)+3 _^ q ^ ^ ^ Q 

jQ T r\{u<ri} J 



IQ T n{u<ri} 

Analogously, it is easy to check that all of the other integrals over {u < i]} on the 
right-hand side of (3.22) are bounded from above by some continuous function, 
h(rf), such that h(rj) — > as r\ — > 0. Therefore, first passing to the limit e — >• 0, and 
afterwards letting 77 — > 0, we easily obtain 



( 6 u n u 2 xxx dxdt < 

1{«>0} 



2- 1 f ( 6 \u x (x,T)\ 2 dx + f 
Jn JQ T n{ 

2 _1 / ( 6 \u 0x \ 2 dx- u n u xxx [2u xx (( 6 ) x + u x (( 6 ) xx }dxdt- 

Jn JQ T n{u>Q} 

[ {vu m u x - Au M u x )(u x ( 6 ) xx dxdt. (3.31) 

JQ T n{u>o} 

Since u(-,t) G C 1 (17) for almost i G [0,T], it is possible to estimate from be- 
low the second term on the left hand side of (3.31) using the generalized Bernis 
inequalities given in Lemma 6.3. As result we obtain: 

2- 1 / ( e \u x (x,T)\ 2 dx + d 7 f C,s{{u n ¥)l + \{u^) xx \ i + (u n ¥) 2 xxx )dxdt + 
Jn Jq t 

do f ( 6 u n u 2 xxx dxdt < 2- 1 [ ( 6 \u 0x (x)\ 2 dx + d 10 [ \( x \ 6 u n+2 dx dt- 

jQ T n{u>0} Jn JQt 

/ u n u xxx [2u xx (C 6 ) x +u x (( 6 ) xx ]dxdt- 

JQ T n{u>a} 

[ {vu m u x - Au M u x )(u x ( 6 ) xx dx dt, (3.32) 

JQ T n{u>Q} 

Next we estimate the terms in the third integral on the right-hand side as in 
[30]; namely, 
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(C )xx dxdt 
rn{u>o} 



/ 

= 6 / {u^u xxx ( 6 ){u^u x Q{u^{b( x +QQ xx ))dxdt 
JQ T n{u>o} 



\Jg T n{«>o} / \JQ T n{«>o} / 

/ u n+2 (5g + (( xx ) 3 dxdt) ' <S [ (u n ul xx + {u^f x )edxdt 

JQt / JQ T n{u>0} 

+ c(S) [ u n+2 (£ + (CCxx) 3 ) dxdt, v<j > 0, 



and 



/ u n u xxx u xx (( 6 ) x dx dt = 6 / (m^w^C 3 ) 

JQ T n{u>o} JQ T n{ti>o} 

"/ n-l 71 — 1 n-4 „ n — 1 ti-4 „\ „"| n n-1 

x 3 u xx H — u 3 u ^ _u 3 u ^J( z j( U 2 3 ( x )dxdt 

1 /2 

\jQ T n{u>o} / L JQ T n(«>o) vn + ^ 



Q T n{«>0} 

1/3 / , \ 1/6 



+ ^(^)V^)V^] (jf u^Cdxdt 



J 

+ c(S) [ u n+2 Cdxdt, V<5 > 0. 

JQt 



<S I (u n u 2 xxx + \(u^) xx \ 3 + (u^f x )C 6 dxdt 

{u>0} 



Using these estimates in (3.32) with 6 = 0*9/6, we obtain (3.21). □ 
4. Finite speed propagation (strong slippage: < n < 2) 

In this section, we consider problem (P) with initial data, Uo, which satisfies 
(2.1) and which also possesses the additional property, 

suppuo C {x < 0}. (4.1) 

Let us introduce the following family of subdomains: 

fi(s) = ft n {x | x > s} Vse(-s,a), Q t (s) = fi(s) x (0,t). (4.2) 

Theorem 5. Let u satisfy (2.1), (4.1), let v = ±1, < A, < n < 2, m < M if 
< A, m<n + 2ifu=l, A = 0, and 

to > if v = —1, to > — i/ i/ = 1, 

and Zef u denote an arbitrary strong nonnegative solution of problem (P), obtained 
as in Theorem 1, which satisfies the local entropy estimate in Theorem 2 or 3. Then 
u possesses the finite speed of propagation property in the sense that there exists a 
continuous function, s(t), such that s(0) = 0, and a positive time T , such that 

suppu(-,t) C Tl\n(s(t)), s(t)<a Vi <T , s(T ) = a. (4.3) 

Remark 4.1. The analysis in the section applies also to (1.4 i )■ 
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Proof. Since by assumption < n < 2, the local entropy estimate, (3.17) or (3.18), 
holds for some positive a < 2 — n. The proof of the finite speed of propagation 
property is based on careful analysis of the properties of solutions satisfying these 
inequalities with positive a. In the stable case, v = —1, such analysis was performed 
in [20] with A — 0. If v = — 1 and A > 0, a similar proof can be given. Therefore 
we restrict our attention here to the unstable case, v = 1. Although we set A = 
for simplicity, all our estimates are also valid for the case v = 1, A > 0. Thus, the 
estimate (3.18) can be written in the form 

, * / ?u 1+a {x,T)dx + Cl f C 4 K + "- 27+1 (« 7 )L + « a+ "^)^^ 
a{a + 1) J n J Qt 

< , 1 - / C 4 «o + " dx + cR, a> 0, (4.4) 
a(a + 1) J n 

where the constant c = max(c 2 , c 3 , c 4 ) and 

R:=R 1 +R 2 +R 3 = f (\C\ + \CC xx \> n+a+1 dxdt+ 

\(( 3 ( x ) x \u a+m+1 dxdt+ [ ( 4 u a+2m - n+1 dxdt. 

Qt JQt 

Note now that due to the constraint m > n/2, the minimal power of u(x, t) in terms 
of the form Q 4 u ai dxdt on the right hand side of (4.4) is strictly greater than 

the minimal power of u in terms of the form J Q ( 4 u a2 dx appearing on the left hand 
side of (4.4); namely, 

a + 2m-n + l>a + l m > —. 

This shall allow us to derive a functional relationship from (4.4), for a suitable set 
of energy functions, which satisfies the conditions of the generalized Stampacchia 
lemma, Lemma 6.2. 

To put (4.4) in an appropriate format, we set ((x) = ( Sy s(x), where ( s ,s(x) is 
the cut-off function defined as 

Csi5 (z) = ^(^i), (4.5) 

where s £ M, 5 > are free parameters, and ip(r) is a nonnegative nondecreasing 
C 2 (R) function such that: 

<p(r) = for r < 0, <p(r) = 1 for r > 1. (4.6) 

We now define three energy functions which are connected with the terms on 
the right hand side of our entropy estimate (4.4), 



J T (s):= u Pl+a+1 dxdt,E T {s) := / u h+a+1 dxdt, 



I t (s) := [ u^ +a+1 dxdt, 



(4.7) 



where (3\ = n, 02 = m, (3 3 = 2m — n. Using (4.4), we shall deduce a system of 
three functional inequalities for Jt(s), Et(s), It(s)- Having set ((x) — ( s ,g(x) in 
(4.4), we obtain after some simple computations 
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dxdt 



sup / u a+1 da; + £>i / [(«^)^ x + (u ! " ± Z ±1 ) 

l(Q.T)Jn(s+S) JQt(s+S) 1 

4 / u a+n+1 dxdt+ 4 / u a+m+1 d:rdi+ 

d JQt(s)\Qt(s+S) ° Jqt(s)\Qt(s+8) 

f 



a+2m—n+l 



'Qt(s) 



+ I ul +a dx. (4.8) 

fi(s) 



Here and throughout the proof, Di denote positive constants which can depend on 
the problem parameters, a, n, to, but not on s, 6, and T. For arbitrary (3 > 0, Z > 
4, s > 0, <5 > 0, such that s + 25 < a, we obtain for a — f3 + a + 1 that 



/ u a dx< f (uCl +s , s ydx= f 



, s +s,s) 7 dx= I vn+S+i dx, (4.9) 

where v — v(x,t) = (u ( l s+s s ) . By the Gagliardo-Nirenberg interpolation 
inequality, 



/ \v x \ A dx 


n + a + l 


f 4( Q + 1) 


Q + l 


JQ(s+8) 




JQ(s+8) 





v™+<*+ l dx < D 3 

) 

where = ^ggg^jy - Combining (4.9) and (4.10) 
u l3+a+1 dx <D A I u a+1 dx 



(4.10) 



/ 



fi(s+2<5) 



n(s+s) 



(l-O)T 
Q + l 



/ l( u " +4 ° +1 )J 4dx + i / 

Jn(s+s) Jn( 



u n+a+1 dx 



Let us suppose that 



6a 



s+<5)\f2(s+2<5) 



< 1. 



(4.11) 



n + a + 1 
or, equivalently, that 

(3 < n + 4(a + 1). 

Then, integrating (4.11) with respect to t and using Holder's inequality, 



(4.12) 



f u P+a+l dxdt < ^yl-^+n gup j /* u Q + 1 (t)dx) 

Jqt(s+28) te(0,T)\Jn(s+8) j 



Qt(s+5)\Q t (s+28) 



i n+a+1 dxdt + j |( U ^ ±i ) a 

JQt(s+S) 



dx dt 



(l-0)«3 + Q + l) 
c + 1 



•« + 4(q + 1) 



< £) rT 1 "+4( Q + l) 



sup 

te(0,T) Jsi(s+(5) 



1 /" 



i n+a+1 dxdt 



+ 
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Qt(s+5)\Qt(s+25) 



(u +4+ ) \ 4 dxdt 



Qt(s+8) 



, (4.13) 



where u = ,f ,., > 0. 

^ n+4(a+l) 

Using the definitions in (4.7) and the a priori estimate (4.8), we deduce from 
(4.13) that 



Ic 



Qt(s+2S) 



Jt(s)- Jt(s + 26) 



S 4 



+ 



E T (s) -E T {s + 6) 



+ It(s)+ f < +1 
Jn(s) 



dx 



(4.14) 



The inequality (4.14) holds for the three values of i = 1, 2, 3, prescribed in (4.7), 
if condition (4.12) holds with (5 = fa, i = 1, 2, 3. These conditions may be written 

as 

l)a + l>0, 2) m< n + 4(a + 1), 3) 2m - n < n + 4(a + 1). 
It is easy to check that all of these conditions are satisfied for some a 6 (a 2 , 2 — n) 
if and only if 

m<6-n. (4.15) 
Thus, if the inequality (4.15) holds, we obtain the following system of functional 
inequalities: 



J T (s + 6) < D 8 T^ 



J T (s) - J T {s + 5) + E T (s) - E T (s + 5) 



S 4 



S 2 

I T (s) + h (s) 



l+Mi 



E T (s + S) < D 9 T 



I T (s + S)< i^ioT— 



Jt(s) - Jt{s + 5) | E T (s) -E T (s + S) 



S 4 



S 2 

I T {s) + h {s) 



1 + M2 



(4.16) 



Jt(s) - Jt{s + 5) | E T (s) - E T (s + 5) 



S 4 



S 2 



1 1+^3 



where ho(s) = J n(s) ul +a dx, and 



4n 



Am 



/'2 = 



I T (s) + h (s) 

4(2to - n) 



n + 4(a + l)' '" A n + 4(a + l)' M3 n + 4(a + l)' 
Due to the boundedness and nonnegativity of u, the following estimates are obvious 



J T (0) < J T := [ u f3l+a+1 dxdt < cT, VT > 0, 
Jq t 

E T (0) <E T := [ u f}2+a+l dxdt < cT, VT > 0, 
Jq t 

J T (0) < I T := / ?/ 3+Q+1 rfsdi < cT, VT > 0, 

'<3t 



(4.17) 



where c is a constant which does not depend on T. The validity of the statement 
of Theorem 5 when inequality (4.15) holds, now follows from (4.16), (4.17), and 
Lemma 6.2, since h (s) = for any s > 0. 

If m > 6 — n, we proceed as follows. We fix m such that § < m < 6 — n. It 
is easy to see that due to the boundedness of the solution u, all of the previous 
estimates in the proof of Theorem 5 remain valid when m is replaced by to. As 
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a result, the system (4.16) is obtained with respect to new energy functions (4.7) 
defined by the values 

Pi =n, 02 = m, 03 = 2 m - m, 

and where 

An 4m 4(2 to — n) 



Ml n + 4(a + l)' M2 n + 4(a + l)' ^ n + 4(a + l)' 

In this manner, the validity of the statement of the theorem for the case to > 6 — n 
again follows from (4.16) and Lemma 6.2. □ 

5. Finite speed propagation (weak slippage: 2 < n < 3) 

In this section we shall again consider problem (P), with initial data, w , which 
satisfies (2.1) as well as the additional property, (4.1). The subdomains, fi(s) and 
Qt(s), are to be understood here to be as defined in (4.2). 

We first prove the following lemma, which provides control on the Lj oc (fl) norm 
of some minimal positive power of the solution under consideration, u(x,t). For 
the sake of simplicity, the results in this section are proven for v = 1 and A = 0, 
though they remain valid for v = — 1 and A > as well. The results here can also 
be readily shown to apply to (1.3) if 2 < n < 3 and to (1.4 i). 

Lemma 5.1. Let v = ±1, < A, 1/2 < n < 3, to > n/2, r\ > ^jp, e > 0, with 
m < M if A > and m<n + 2ifu = l, A = 0. Then there exists a positive 
constant c, depending on n, m, 77, e only, such that any nonnegative strong solution 
u of problem (P) satisfies 

( u{x,Ty +1 C, A dx<e( [ ( 6 u n u 2 xxx dxdt+ f \( x \ 6 u n+2 dx dt) + 

Jn \JQ T n{u>o} Jq t J 

i n+2r >\c x \ 2 + u ^r 1 - ? + u ™+"+ i icc,i 2 ] dxdtj + 

f ul +1 C, 4 dx + c f vT+^-^dxdt, (5.1) 

JQ, ^QtHsuppC 

for arbitrary nonnegative ( G C 2 ([—a,a\). 

Proof. The proof here follows that of Lemma 5.2 in [25]. Let us set <p — cp in the 
integral identity in (2.2), where ip is the test function 

v = -hCHu + 1 )\ 7 >o, 

where {lg(t)} C C^°(0,T) and Is — > X(o,t) as <5 — > 0. After some simple computa- 
tions, we obtain 



7 K 



j 



(is)tc 4 (m+ 7 ) 1 t ' +1 dxdt^ / u m u ^(( u + 7 )"c 4 ) :c ^^+ 



/ ls(( 4 ) x u n (u + j) n u xxx dxdt+ 



"I 

JQ T n{u 



Q T n{u>o} 

l s ^u n (u + j) r) - 1 u x u xxx dxdt := A 1 +A 2 + A 3 . (5.2) 



>o} 
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The terms A2, A3 can be estimated as in [25, Lemma 5.2]. For any e > 0, 

\A 2 \<e[ (?u n ul xx dxdt + Ci{e) [ \( x \ 2 u n (u + 7 ) 2 " dx dt, 

JQ T n{u>0} JQt 

\A 3 \<e( [ ( 6 u n u 2 xxx dxdt+ [ \( x \ 6 u n+2 dxdt) + 

\JQ T n{u>0} JQt J 

C 2 (e) f (u + 1 ) n+ ^- 1 dxdt. 

J QrHsupp C, 

Here Cj denote constants which may depend on m, n, r\, and on e if indicated, but 
which are independent of 7 and 8. 
Let us now estimate A\. 

Ai= f r l u m {u + 1 ) n - 1 u 2 x C, A l s dxdt+ 

JQt 

f 4u m (u + >y) r iu x C 3 {xl8dxdt:=A? ) + A™ . (5.3) 

JQt 

Since m — +77 — 1 > 0, it follows from Young's inequality and Lemma 6.4 that 
for any e > 0, 

\A^\<e[ u 6 x u n - 4 ( 6 dxdt + C 3 (e) [ {u + 7 )i( m -^+"- 1 )c 3 dx dt < 

JQ T n{u>0} JQt 

( 6 u n u 2 xxx dxdt + C 4 f £u n+2 dxdt) + 

l T n{u>0} JQt J 

C 3 (e) / (u + j) 3m+3 l +1 ~" S 3 dxdt. (5.4) 



(2) 

With regard to , we have by Young's inequality 

|4 2) | < \A™\ + C 6 [ {u + 1 ) m+r i +1 \tt x \ 2 dxdt. (5.5) 

JQt 

Thus all the integrals in (5.2) are uniformly bounded with respect to the pa- 
rameters 8 > 0, 7 > 0. Therefore, collecting the estimates obtained for the terms 
Ai,i — 1,...,3, passing to the limit 8 —¥ 0, and then to the limit 7—^0, the 
estimate (5.1) follows. □ 

Theorem 6. Let u satisfy (2.1), (4.1), let v = ±1, < A, 1/2 < n < 3, m > n/2, 
m < M if A > 0, and m<n + 2ifv=l, A — 0, and let u denote an arbitrary 
strong nonnegative solution of problem (P) obtained as in Theorem 1, which satisfies 
the local entropy estimate in Theorem 2 or 3. Then u possesses the finite speed of 
propagation property in the sense of Theorem 5. 

Proof. Let us consider the local energy estimate (3.21) obtained in Theorem 4, and 
estimate the third term on the right hand side, setting v = 1 for simplicity, 

B := / u m u x (u x ( 6 ) xx dxdt = 

JQ T n{u>o} 

/ u m u x u xxx ( 6 dx dt + 12 / u m u x u xx ( 5 ( x dx dt+ 

JQ T n{u>o} JQ T n{u>o} 
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Jq t 

Using Young's inequality, we obtain that for any e > 



6 / u ro u2(C 8 C*)x dxdt := B 1 + B 2 + B 3 . 

n{u>o} 



\Bi\< 



n (ti-4) „ n-4 a 

u 2 u xxx u x u e u m 2 6 dx dt 

Q T n{u>Q} 



< 



■ f (u n U 2 xxx + <W"- 4 )C 6 dx dt + D 1 (6) f u 3 m -2n+2 c 6 

jQ T nfu>0} JQt 



dxdt. 



lQ T n{u>o} 

Here and in the sequel Di denote constants which may depend on m, n, rj, and on 
e if indicated, but which are independent of 5 and s. Similarly, we may estimate 

B 2 <e[ (u^u"- 1 +ulu n - 4 X 6 dxdt + D 2 (e) [ u 2m - n+2 \( 2 ( x \ 2 dxdt, 

JQ T ri{u>0} JQt 

B 3 <e[ u 6 x u n - 4 C 6 dxdt + D 3 (e) [ \(( 3 ( x ) x \i dxdt. 

JQ T n{u>0} JQt 

Using these estimates in (3.21), we obtain that for e > sufficiently small 
/ \ Ux (x,T)\ 2 ?dx+^ [ C%u^t + (u^)l x + {u^) 2 xxx ]dxdt+ 

Jn z JQt 



dw f 

2 JQt 



( 6 u n u 2 xxx dxdt < 
n{u>o} 



f \u 0x \ 2 C 6 dx + d 9 [ (\( x \ 6 + \tt xx \ 3 )\u\ n+2 dxdt+ 

Jn JQt 

dJ [u 3m - 2 "+ 2 C 6 + u 2m -" +2 |C 2 C,l 2 + «'"^l(C 3 C,),| i ]^^. (5.6) 

JQt 

Assuming that r\ > i^p, summing the inequalities (5.6), (5.1), and taking e > to 
be sufficiently small, we obtain: 

/ \u(x,T)\^(Ux+ f \ Ux (x,T)\ 2 ( (i dx + ^ f C 6 \(u^) xxx \ 2 dxdt 
Jn Jn 4 JQ T n{u>o) 

< / \u 0x \ 2 ( 6 dx + / |w r +1 C 4 dx + D 5 R, 
Jn Jn 

R:= [ u 3m - 2n+2 ( e dxdt+ f u m+3 2 +1 ~ n ( 3 dxdt+ 
JQt JQt 

j u n+ ^d x dt+( (|Cx| 6 + |CCzxl 3 K +2 <fa^+ 

•^QtHsuppc; JQt 

f u 2m - n+2 \( 2 ( x \ 2 dxdt+ f u n+2r >\( x \ 2 dxdt+ 

JQt JQt 

f u m +^ +1 \a x \ 2 dxdt+ f u^f^^C^dxdt. (5.7) 

JQt JQt 

Let us take C,{x) as the function ( Sy $ from (4.5). It then follows from (5.7) that 

sup / \u(x,t)\ 7 > +1 dx + sup / \u x (x, t)\ 2 dx 

te(o,T)Jn(s+s) te(o,T)Jn(s+s) 
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dw 
4 



R:= 



/ \(u S ¥-) xxx \ 2 dxdt < (\u 0x \ 2 dx+ \u \ v+1 )dx + D 6 R, 

Jq t (s+s) Jn(s) 

f u 3m - 3n+2 dx dt + [ M 3 "+ 3 g+ 1 -" da . dt+ /* u ™+ 3 ')- 1 ( ix ( it+ 

JQt{s) ->Qt(s) JQt(s) 

<T 6 f u n+2 dxdt+ 
5- 2 [ u 2m - n+2 dx dt + S- 2 [ u n+2r < dxdt + S- 2 [ u rn+r ' +1 dxdt+ 

JQt{s) JQt{s) JQt(s) 

S~ 3 f u^^ 1 dxdt:=J2 S ~ Xl [ u^dxdt. (5.8) 

JQt{s) i=1 JQt(s) 

We shall use the inequality (5.8) to derive a system of functional relationships, 
which will allow us to implement the Stampacchia Lemma, Lemma 6.2. In order to 
undertake similar computations to those undertaken in §4, it suffices to guarantee 
that the following inequalities hold: 

&>l + m i = l,2,...,8. (5.9) 

First we ensure that 

n . . 

e 3 = n + 3ry-l>l + r7^?7>l--:= ry mm . (5.10) 

Next we deduce a restriction on m by considering 

£i = 3ra-2n + 2>l + r?<^ -(2m - ra) + 1 + (l - -J > 1 + 77. (5.11) 

Together, (5.10) and (5.11) yield that 

n ■ =1-- 
There exists 77 satisfying (5.12) iff 



n 3 



2m-n>0. (5.13) 

Next it is easy to see that 

. 3m + 3r?+l-n (n + 3ry - 1) + (3m - 2n + 2) 
6 = 2 = 2 • 

Therefore the inequality 

£2 > l + r? (5.14) 

follows from (5.10) and (5.11). It is easy to check that the other inequalities in 
(5.9) can be satisfied by an appropriate choice of 77 if conditions (5.10) and (5.11) 
are satisfied. 

As result of (5.9), the following Gagliardo-Nirenberg interpolation inequalities 
hold for i = 1,2,..., 8, 

<>iU 

dx 



f u^dx<D 7 l[ |((ci»' 

J(l(s+28) \jQ(s+S) 

x (/ &» n+ld A " . (5-15) 

\JQ(s+5) I 



T) + l 
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where &i = ^^ n ^ , i — 1, 2, . . . , 8. We now wish to guarantee that 

^<1, i = l,2,...,8. (5.16) 
From the definition of 9 i} it follows that (5.16) holds iff 

V> €i ~"~ 8 , i = l,2,...,8. (5.17) 
6 

It is easy to check that all the inequalities in (5.9) and (5.17) hold for some 77 in 
the interval (5.12) if condition (5.13) is satisfied. 

Therefore we deduce from (5.8), (5.15) the following inequalities 



f dxdt < DsT 1 -^ (y25~ x > [ u^dxdt+Ho 

JQt(s+8) V =1 Jqt(s) 



1+ 

(*)' 



5>7+n+7 



(5.18) 



where H Q (s) := J n(s) (wo + ' 7 + \u 0x \ 2 ) dx = 0, for all < s < s + S < a. From (5.18) 
and Lemma 6.2, the conclusion of Theorem 6 now follows. □ 

Remark 5.2. The method used here for the qualitative analysis of conditions which 
guarantee the finite speed of propagation property, should be possible to modify, mod- 
ulo various technical difficulties, to yield quantitative predictions for propagation 
rates, [25, 30, 45]. We hope to perform such an analysis in a forthcoming paper. 

6. Appendix A 

Lemma 6.1. (Stampacchia's Lemma [48] Assume that a given non-negative, 
nonincreasing function g : [0, 00) — ► R\_ satisfies 

g(s + S)^co(S- a g( S )Y 3 , Vs > 0, 5 > 0, 

for some real numbers Co > 0, a > 0, [3 > 1. Then 

g(s )=0, s :=2^T( Co 5 (0)' 3 - 1 )^. 

For given (f3\, ■ ■ ■ , (3k) £ we denote 

k k 

Lemma 6.2. (A Stampacchia Lemma for Systems [25]). Assume that k given 
nonnegative non-increasing functions, g; t : [0, 00) — > [0, 00), i = 1, • • ■ , k, satisfy 

/ k \ A 
gi(s + S) \V,5- ai 9i (s)J , Vs>0, S >0, i = !,-■■ ,k, 

for some real numbers Ci > 0, > 1; di ^ 0, i = 1, • • • , k, and a.i > for 
i = 1, • • • , I, I > 1. Denote 



i=l 
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and assume that the function Q(s) defined by 



Q(s) 



i=i+i 

o, 



k - / — \ 



i/ Z < fc, 



i/ Z = fc, 

satisfies the constraint Q(si) < 1 at some point S\ ^ 0. TZien ZZiere exists a positive 
constant C > 1, which depends on k,l,Oi,j3i, and Q(s\), such that 

9i(so) = 0, V i = 1, ••■ ,Z, 

w/iere 

*o ^.i+cj;^^) 1 -*^))*- 1 

Remark 6.3. Ifl = k, then s\ = in Lemma 6.2 

Lemma 6.4. (Generalized Bernis inequalities for periodic functions [25]/ 
Let Q, = (a, b) C M 1 , n G (|, 3), and Zet u G C^fj) n #f oc ({u > 0}) &e an arbitrary 
|0| periodic function such that 



u ^ 0, 



u x e L 2 {Q), 



' dx < oo. 



{u>0} 



Tften bx g W 1 ' 6 ^), G iy 2 < 3 (f7) ; G T^ 3 ' 2 (fi) := H 3 (0), and there 

exists a positive constant, C > 1, which depends on n only, such that 



dx 



/ <•(!(.*) r + |(-*) r + |(-*) 

y\ / X V /XX V / Xi 

J{«>0} ./n 
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